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Abstract: - For the semilinear elliptic equation A u + p(x)u~' 1 = 0, x £ R n , 
n > 3, 7 > 0, we show via the barrier method the existence of a positive 
entire solution behaving like |a,*| 2-n near oo. 

1 Introduction 

We study the singular semilinear elliptic equation 
(1) Au + p(x)u~' 1 = 0 

in R n . This type of equation arises in the boundary layer theory of viscous 
fluids [3,4]. From the results of Fulks and Maybee [8], Crandall, Rabinowith, 
and Tartar [5], Gomes [9], and recently Lazer and McKenna [14], it follows 
that (1) has a unique classical solution within a bounded domain Q, where 
p(x) is a sufficiently regular function which is positive on Cl. Kusano and 
Swanson [12] gave the existence proof on exterior domains. As for the ex- 
istence of entire solutions, not much is known. Edelson [7], Kusano and 
Swanson [13] have been able to show the existence of entire solutions of (1) 
with 7 G (0, 1), and p(x) sufficiently regular. In this paper we show via the 
upper and lower solution method, which is also referred to as the barrier 
method, that (1) has a bounded positive entire solution vanishing at oo in 
R n for n > 3 and all 7 > 0. 

The author learned after this paper was finished that a similar result was 
given earlier by R. Dalmasso [6], but by a different approach. 
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2 Preliminaries 



We first state the theorem by Kusano and Swanson [ 13 ] for the case 0 < 7 < 
1. 

Lemma 1. Equation (1) has an entire bounded positive solution u(x) in R n 
for n > 3 ,and |a:| n_2 u(a;) is bounded and bounded away from zero near 00 if 
p(x) satisfies the following conditions: 

(HI) p{x) G C? oc (R n ), n > 3 , p(x) > 0 , x G R n \ { 0 }, 

(H2) 3 C > 0, such that C<f>(\x\) < p(x) < <f>(\x\), <j>\ x) = max\ x \ =t p(x), 0 < 
t < 00, 

(HS) f^° t n - 1+ ^ n -^<f>(t)dt < 00. 

The term ’’entire” has often been used for solutions of equation (1) in 
R n . To avoid confusion with the traditional definition for entire functions, 
we use the term ”C 2+a -e ntire”. A C 2+a -entire solution of (1) is defined to 
be a function u(x) G Ci+ a (R n ) that satisfies ( 1 ) pointwise in R n . 

The method that we shall be using heavily in our proof is the so-called 
barrier method, or upper-lower solution method. 

We consider the elliptic boundary value problem 

. . f Lu + f(x,u) = 0 in D 

| Bu = adu/du + bu = g on dD 

where D is a smoothly bounded domain in R N and v = (17, • • • , i/ n ) is a 
smoothly varying outward normal vector field on dD which is of class C 2+a , 
while a and b are positive constants. We also assume that / G C a and that 
g has an extension g to the interior of D such that g G C 2+a . 

An upper solution to the above problem is a function <^> satisfying 

f L<f> + f(x, (f>) < 0 in D 
\ Bf> > g on dD. 

A lower solution to the above problem is a function xf> satisfying 

f Li[> + f(x,xj>) >0 in D 
Btf < g on dD. 

We assume that dD,f,g, and the coefficients of L are smooth in what 
follows. 
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Lemma 2 . (Theorem 2.3.1 of [16]) Let <p be an upper solution and ip a 
lower solution with ip < <p on D . Then there exists a solution u to the above 
boundary value problem with ip < u < (j). 

We consider the following example: 

j u " ' + Xu — u 3 = 0 x € (0, 7r) 

| u = 0 x = 0,ir 

By the above theorem, if A > 1 , then the problem has at least three solutions. 

Actually, u = esinrr with e small is a lower solution, and u = Rx '/ 2 with 
R large is an upper solution. Therefore there exists a solution u such that 
u < u < u in ( 0 , 7r) . Clearly — u and 0 are also solutions to this problem. 

The following lemma on the barrier method for D = R n is due to Ni 
[15] in 1982. A special case was proved earlier by Ako and Kusano [ 1 ] in 
1964. The proof is standard. Using the well known result on the upper-lower 
solution approach in bounded regions (see Sattinger [16]), we first solve the 
equation 

Lu + F(x,u) = 0 

on Br. Then by letting, R — * oo, we obtain a solution on R n by a diagonal 
process. 

Lemma 3. Let uj > u 2 in R n be such that 

/ox f Lui + f(x,Ui) < 0 

^ ' \ Lu 2 + f{x,u 2 ) > 0 

where f is locally Holder continuous in ( x,u ) and locally Lipschitz in u, and 
L is an elliptic operator of second order. Then there exists a solution u of 
Lu + f(x, u) = 0 with Ui > u > u 2 . 

3 Main Result 

Theorem 1 . Under the same conditions as given in Lemma 1, the equation 
(1) has a C 2+a -entire positive solution in R N ,N > 3, vanishing at oo at the 
rate of at least |ar with some q € ( 0 , 1 ) for any 7 > 0 . 
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The difficulty in constructing the proof is to find an appropriate upper 
solution for equation (1). In order to use the barrier method we first study 
the nonsingular equation 



A u + p{x )\8 + u] 7 = 0. 



For each fixed 7 there corresponds a solution u y (x). Letting 7 — > 00, we 
show that the limiting function is the desired solution. 

Proof: By Lemma 1, for 7 = 71 £ (0,1), equation (1) has a C 2+0, -entire 
positive solution in R n , n > 3, vanishing at 00 at the rate r 2-n . We 

claim that u = cu q is an upper solution of the equation (1) for 7 > 1, where 



- 1 T 7i 1 

9 < ~rr~ < 

1 + 7 

^fl+71 -?(!+')') _j_ 

c > ( ) 1 +t , M = max x€ Rn\u(x)\. 



In fact: 



A5+d£) 

IP 

cq(q — l)u ?-2 |Vu| 2 — cqu q ~ l p(x)u ~ J1 + p{x)c ~ 1 u~ iq 
< —cqu q ~ l p(x)u ~~ 11 + p(x)c~' 1 u~~ ,q 



p(x) M C 1+7 <? 






( 1 - 



u 



1+71 -?(!+■>) 



) 



< 2 ^ 1(1 






1+71-9(1+7) 



U 



) 



<£W (1 - 1) = 0. 



Let 8 be a fixed positive number. We then observe that u is an upper 
solution of the equation 



(4) A u(x) + p(x)[u(x) T 5] 1 = 0, x G R n . 

u = 0 is a lower solution of (4). Since u = cu\ > 0, u > u in R n . By Lemma 
2, (4) has a solution u such that u<u<u. 
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For 8 < 8, u is a lower solution of (4) with 8 = 8. Lemma 2 then implies 
that (4) has a solution u for 8 = 8 such that u < u < u. 

Let {^n } i° be a sequence of strictly decreasing positive numbers, and let 
u n (x) be a smooth positive solution of (4) when 8 = 8 n . From the construc- 
tion of our lower solutions, it is clear that u n (x) > u„_i(x) for all n. So 
limn—oc u n (x ) = u(x) exists for all x € R n and 

(5) u < u < u 



for x (E R n ■ 

We can now assert that u 6 C 2+a (R n ) and that 

(6) An + p(x)u~ 1 = 0 

for x € R n - This follows from more or less standard arguments. 

Let x 0 £ R n and r > 0. We consider the ball of radius r cent ered at x 0 , 
B(x 0 ,r) in R n . Let T be a C°° function which is equal to 1 on B(x 0 ,r/2 ) 
and equal to 0 off B(x 0 ,r). We have 

A ($u n ) = 2VT • Vu n + p n 

for n > 1, where p n is a term whose L°° norm is bounded independently of 
n. Therefore for n > 1 we have 

T?i n A(Tu n ) = b nj-\ + 9n, 

j = 1 V X 3 

where b n j,j = l,---,n and q n are terms bounded independently of n for 
n > 1. Integrating the above equation, we have that there exist constants 
Cj > 0 and c-i > 0 independent of n such that 

/ |Vu n | 2 dx < c x { f |Vu n | 2 dx) 1//2 + c 2 . 

Jb(xo,t) Jb(x 0 ,t) 

From this, it follows that the L 2 (B(x 0 , r))-norm of IV'Lunl is bounded inde- 
pendently of n. Hence, the L 2 (B(x 0 , r/2))-norm of |Vu n | is bou nded inde- 
pendently of n. Let Tj be a C°° function which is equal to 1 on B(x a ,r/ 4) 
and equal to 0 iff B(x 0 ,r /2). We have for n > 1, 



— 2VT] • Vu n + Pi„, 
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where p ln is a term whose L°°(B(x 0 ,r / 2))- norm is bounded independently 
of n. From standard elliptic theory, the W 2,2 (B(x 0 ,r / 2))-norm of is 

also bounded independently of n and hence, the W 2,2 (B(x 0 , r/4))-norm of 
u n is bounded independently of n. Since the W 1,2 (B(x 0 ,r / 4))-norm of the 
components of Vu n are bounded independently of n, it follows from the 
Sobolev embedding theorem that if q = 2n/(n — 2) > 2 for n > 2 and in 
addition if q > 2 is arbitrary for n < 2, then the L q (B(x 0 , r/4))-norm of |it n | 
is bounded in dependently of n. Let $2 be a C°° function which is equal to 1 
on B(x 0 ,r/ 8) and equal to 0 iff J5(x e ,r/4). We have for n > 1, 

A($ 2 u n ) = 2Vtf 2 • Vu n + p 2n , 

where p 2n is a term whose L°°(B(x 0 ,r/ 4))- norm is bounded independently of 
n. Since the right hand side of the above equation is bounded in L q (B(x 0 ,r / 4)) 
independently of n, the W 2,q (B(x 0: r/4))-norm of ^ 2 u n is also bounded inde- 
pendently of n. Hence, the VF 2,, (R(:r 0 ,r/8))-norm of u n is bounded indepen- 
dently of n. Continuing this line of reasoning, after a finite number of steps, 
we find a number rq > 0 and <jq > n/( 1 — a) such that the W 2,qi (B(x 0 ,ri))- 
norm of u n is bounded independently of n. Hence, there is a subsequence 
of {u n }J°, which we may assume is the sequence itself, which converges in 
C 1+a (B(x 0 ,ri)). If 0 is a C°° function which is equal to 1 on (B(x 0 ,ri/ 2)) 
and 0 off B(x 0 ,r 1 ), then 

A (6u n ) = 2V0 • Vu n + p n , where p n = 6Au n + u n A6. 

The right-hand side of the above equation conve rges in C°^( B(x 0 ,ri)). Hence 
by Schauder theo ry, {fa n }” co nverges in C 2+a (B(x 0 , ri)) and thus {u n }i° 
converges in C 2+a (B(x 0 ,ri/2)). Since x a was arbitrary, this shows that u 6 
C 2+a (R n ). Clearly (6) holds. 

4 Some remarks 
Remark 1: 

For n = 1, the properties of positive solutions of equation (1) have been 
studied by Taliaferro [17], and Gatica [10]. For n = 2, no entire positive 
solution of equation (1) exists regardless of its asymptotic behavior at oo 
(see [13]). 

Remark 2: 
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It is observed by Callegari, Friedman and Nachman [2], [3,4] that if the 
partial differential equations describing the boundary layer behind a rarefac- 
tion or shock wave (with viscosity proportional to the temperature) traveling 
down, and perpendicular to, a flat plate are written in terms of a stream func- 
tion and a similarity variable the following Blasius-type equation emerges 
[18]. 

nv) + mnv) = o, 

where 

m = 0. /'( 0) = K, /'( oo) = 1. 

Here, 0 < K < 1, for rarefaction waves and, 1 < K < 6, for shock waves. 
(I( = 0 corresponds to the classical Blasius problem.) Adopting the Crocco 
variables 

* = /'(*?), g = f"(v) 

results in the system 

gg" + x = 0, 
g\I<) = 0, g( 1) = 0, 

which falls into the class of equation discussed in this paper. 
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